Spin dynamics and magnetic correlation length in two-dimensional quantum 

Heisenberg antiferromagnets 
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The correlated spin dynamics and the temperature dependence of the correlation length £(T) in 
two-dimensional quantum (S = 1/2) Heisenberg antiferromagnets (2DQHAF) on square lattice are 
discussed in the light of experimental results of proton spin lattice relaxation in copper formiate 
tetradeuterate (CFTD). In this compound the exchange constant is much smaller than the one in 
recently studied 2DQHAF, such as La2Cu04 and Sr2Cu02Ci2. Thus the spin dynamics can be 
probed in detail over a wider temperature range. The NMR relaxation rates turn out in excellent 
agreement with a theoretical mode-coupling calculation. The deduced temperature behavior of £(T) 
is in agreement with high-temperature expansions, quantum Monte Carlo simulations and the pure 
quantum self-consistent harmonic approximation. Contrary to the predictions of the theories based 
on the Non-Linear a Model, no evidence of crossover between different quantum regimes is observed. 



In recent years, it has been argued that the pairing 
mechanism leading to high-temperature superconductiv- 
ity can occur in the presence of quantum critical fluctua- 
tions, as the ones expected in two-dimensional quantum 
(S = 1/2) Heisenberg antiferromagnets (2DQHAF) |],§. 
This early observation has triggered a novel interest on 
the correlated spin dynamics in 2DQHAF. Several stud- 
ies have been carried out on the behavior of the magnetic 
correlation length £ as a function of temperature and/or 
disorder, and on the possibility of driving quantum phase 
transitions by hole doping. In particular, for a given de- 
gree of disorder (i.e. doping) the temperature is expected 
to cause the transition to a regime with peculiar quan- 
tum spin excitations, which could be responsible for novel 
pairing mechanisms j^]. In spite of an enormous amount 
of experimental and theoretical efforts this issue is still 
controversial |^|||. Thus, it is believed that it is necessary 
first to unravel the problem of the spin dynamics and of 
the temperature dependence of £(T) in pure 2DQHAF. 

The 2DQHAF with nearest-neighbor interaction on 
square lattice, is described by the Hamiltonian 
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Sj+a , 



(1) 



with J > and a = (0, ±a), (±a, 0), a being the lattice 
constant. The 2DQHAF has been commonly described 
in terms of the 2-dimensional Quantum Non-Linear a 
Model (QNLcrM) f|, with the action given by 
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dxdr (|Vn| 2 -|9 r n| 2 ); |n| 2 = 1 , (2) 



where n(x) is a unitary 3D vector field, g = c s A/p s and 
u = c s A/T are the coupling and the imaginary-time cut- 
off respectively. This mapping is rigorous only in the 



large-spin limit, where p s and c s are the spin stiffness 
and the spin- wave velocity of the correspondent magnetic 
model. This approach can be extended to S = 1/2 pro- 
vided that c s and p s are introduced as fitting parameters. 
Then a rich and interesting phase diagram can be devised 
Accordingly, the QNLcrM has arisen to a great 
popularity, both for basic quantum magnetism studies as 
well as a good candidate to describe the modifications 
in the microscopic properties of a 2DQHAF when it be- 
comes a high-T c superconductor upon charge doping. 

One of the most striking predictions is the presence 
of the so called quantum criticality, in a well defined 
temperature window. Below a finite value of g, by in- 
creasing the temperature, the system should cross from 
the renormalized classical to the quantum critical regime 
and eventually to the classical behavior |^.|[ . 

It was suggested that this scenario could be probed by 
studying the temperature behavior of the static correla- 
tion length £(T). The latter can be directly estimated 
by neutron scattering or by the longitudinal NMR 

relaxation rate, which provides an indirect measure of 
f @,H. While for small £ (~ a), neutron scattering 
studies become difficult and accuracy problems arise, it 
has been proved Jl^] that NMR relaxation, driven by the 
correlated spin dynamics, can be a suitable tool to derive 
the temperature and doping dependence of the magnetic 
correlation length. 

The behavior of £, according to the QNLcM approach, 
considerably changes in the various regimes ■ Until 

now, most of the discussion on the correlation length in 
2DQHAF with 5 = 1/2 has dealt with La 2 Cu0 4 and 
S^CuC^C^. However, in these compounds the large 
value of J ~ 1500 K prevents one to obtain experimental 
results above T / J ~ 0.6 and thus a detailed comparison 
with the theoretical predictions is difficult. As a mat- 
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ter of fact the region where the quantum critical regime 
is expected, according to the QNLerM approach, should 
be 0.4 < T/J < 0.8. Neutron scattering experiments in 
antiferromagnets with larger S have been recently car- 
ried out with good accuracy Jl(],[ll| and no evidence of 
crossover has been found in a large range of T/J. In- 
deed these data are in agreement with high-temperature 
expansion (HTE) results Jiij , quantum Monte Carlo 
(QMC) data |l5| and with the results of the pure quan- 
tum self-consistent harmonic approximation (PQSCHA) 
recently obtained by some of us The latter ap- 

proach is reliable for any realistic temperature range for 
S > 1, and it gives good results starting from T/J > 0.35 
for S = 1/2 0. 

The advantage of CFTD is that besides being a good 
prototype of S = 1/2 2DQHAF on the square lattice 
(the intcrlaycr interaction constant is less than 10~ 4 J), 
still J is rather small (J ~ 80 K) so that the experi- 
ments can be carried out in a wide range of T/J. Thus 
the classical behavior is approached and a comprehen- 
sive analysis becomes possible. In particular, one can 
explore also the temperature range where the correlation 
length is of the same order of the lattice constants a. 1 H 
NMR experiments have been performed on a 5 x 5 x 4 
mm 3 single crystal of copper formiate tetradeuterate 
(Cu(HC0 2 )2-4D 2 0) Q. The spin-lattice relaxation rate 
1/Ti has been measured in a 8 kGauss magnetic field 
(<?/isii <C J) applied perpendicular to the ab plane, by 
using standard pulse sequences. In Fig.l the experimen- 
tal data for 1/Ti are reported. 

The relaxation rates can be calculated, through the 
mode-coupling (MC) approach |T^-p>l[. First, we calcu- 
late the normalized Kubo relaxation function 63l : 



1/Ti = ^J2S k (cj)[A 2 + B 2 (cos(k x a) + cos(k y b))], 
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fik(*) = I (S£(t + i\)S2)d\, 
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by solving the following integro-differential equa- 
tion 
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where Jk = J^ a exp(k-a), and A is a temperature- 
dependent parameter which fixes the static susceptibility 
of the system and can be directly related to £ (T) . 

By means of the spectral theorem the dynamic cross 
section Sk(u>) can be obtained from the Fourier transform 
of -Fk(t). Thus the longitudinal relaxation time is given 
by pi: 



(5) 



where 7 = 2ir ■ 42.576 • 10 2 s _1 Gauss _1 . The form fac- 
tors [A 2 + B 2 (cos(k x a) + cos(k y b))] have been derived 
by assuming the hyperfine coupling of 1 H nuclei with 
the two Cu 2+ nearest neighbors only. The coupling con- 
stants have been estimated from the rotation pattern by 
rotating the crystal around the a- axis in a field of 94 
kGauss. In this way the transferred hyperfine interac- 
tion constant was directly obtained. The dominant dipo- 
lar part was determined through lattice sums (details will 
be published elsewhere p4[). The two coupling constants 
turned out A = 2.31 kGauss and B = 1.46 kGauss. 

The analysis in terms of the theoretical models requires 
the knowledge of the exchange constant J. Ad-hoc neu- 
tron and light scattering experiments, leading to the best 
estimate of J , have been performed in La2Cu04 p5| . 
Analogous methods to evaluate J in CFTD are under 
way j2f|. At the moment in CFTD we can assume a 
value around J = 80 if. Estimates based on the veloc- 
ity of the spin waves f27j] are in the same range, as well 
as old data |2§] (within 10%), in agreement with the dis- 
persion curve determined by neutron scattering near the 
zone-boundary |2lf |. 

The MC scheme requires the static quantities as an 
external input. For this purpose the simplest and self- 
consistent approach is to use those pertaining to the two- 
dimensional spherical model |2lJ . An alternative method 
is to rely on the results obtained for £(T) by the PQSCHA 
to determine the parameter A appearing in Eq. (^), for 
T/J> 0.35. The two theoretical curves are reported in 
Fig.l. Furthermore we have used the more precise QMC 
data for £(T) 15 to calculate the relaxation rate at lower 
temperatures, down to the value of T ~ 20 K, where 
three dimensional ordering effects are revealed by the 
NMR experimental findings. The relaxation rate 1/Ti 
should be weakly temperature dependent in the quan- 
tum critical regime || , without any observable minimum 
around T/J = 0.5 (T = 40 K). Instead the experimen- 
tal data do present such a minimum, which is completely 
reproduced by our approach. Moreover, in a wide tem- 
perature range the behavior of £(T) is also consistent 
with the renormalized classical regime [^9| . The mini- 
mum in 1/Ti is related to the k dependence of the hy- 
perfine form factor in Eq. (^J), which is strongly peaked 
at the center of the BZ but non zero at the AF wave- 
vector k = [%/ a,ir/a). One has to notice that in the 
temperature range far from the transition the spherical 
model seems to be a good approximation, as shown also 
by the temperature-dependence of the correlation length 
in Fig.2. 

The remarkable agreement between MC calculations 
and the experimental results, allows us to derive the cor- 
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relation length £ (T) from the experimental data through 
an inversion of the mode-coupling approach. The results 
are shown in Fig. 2 together with the theoretical predic- 
tions. The success of the mode-coupling approach gives 
further support to the validity of the scaling hypothesis, 
and the correlation length can be extracted with satis- 
factory accuracy from the relaxation rate, by means of a 
procedure which has already yielded reliable quantitative 
estimates in La2Cu04 (see Refs. Jl^ , p^ | for details). In- 
deed, when £(T) ^S> a classical scaling arguments for the 
generalized susceptibility and the decay constants lead to 
the following equation Jl2[: 
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thus establishing a one-to-one correspondence between 
1/Ti and £(T). In Eq. (||) the 2-d wavevector q starts 
from the BZ boundary (ir/a, 7r/a), corresponding to the 
AF ordering wavevector; Lu e is the Heisenberg exchange 
frequency describing the fluctuations in the limit of infi- 
nite temperature, z = 1 the dynamical scaling exponent, 
e = 0.33 takes into account the reduction of the ampli- 
tude due to quantum fluctuations 1 30 1 and (3 is a normal- 
ization factor preserving the sum rule for the amplitude 
of the collective modes. Also the data for £(T) obtained 
by means of this procedure are reported in Fig. 2. It is 
apparent that the values for £(T) are close to the ones de- 
duced by mode-coupling approach, the main differences 
being at low values of £, as expected. 

From the comparison with the theoretical estimates, 
the permanence of the system in the renormalized classi- 
cal regime, on moving towards the classical one at highest 
temperatures, is shown. No evidence of crossover to the 
quantum critical regime appears, as is clearly shown in 
the inset of Fig. 2. Finally we point out that although we 
used a value of J = 80 K in all our calculations, the main 
conclusions derived regarding the absence of a crossover 
towards a quantum critical regime do not depend on J. 

Summarizing, we have used a single crystal of a proper 
2DQHAF which allows one to analyze the behavior of the 
in-plane correlation length £ (T) over a wide temperature 
range in the light of several theoretical models. £(T) has 
been derived from 1 H NMR relaxation rates, having ver- 
ified also that the mode-coupling theory is able to repro- 
duce successfully the temperature dependence of 1/Ti. It 
turns out that £(T) is well described by QMC, PQSCHA 
and HTE. In particular, at variance with other recent 
theories for 2DQHAF, no crossover to the quantum crit- 
ical regime occurs, up to temperatures T/J ~ 1.5, where 
the classical behavior has already set on. Finally, Eq. 6, 
which is based on the validity of the scaling hypothesis, is 
shown to yield reliable estimates of 1/Ti even for £ ~ a. 
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FIG. 1. Proton spin-lattice relaxation 1/Ti as a function of 
temperature. Squares: experimental data; dashed-line: theo- 
retical mode-coupling (MC) results using the static quantities 
of the spherical model; full line: MC results using the correla- 
tion length given by the PQSCHA for T/J > 0.35; filled cir- 
cles: MC results using the correlation length given by QMC 
simulations for T/J < 0.35 [§. 

FIG. 2. Correlation length £ vs temperature. Squares: ex- 
perimental data deduced from 1/Ti by inversion of Eq. (|^); 
circles: experimental data deduced from 1/Ti by inver- 
sion of mode-coupling results leaving A as a free parameter; 
dashed-line: spherical model; full line: PQSCHA; triangles: 
QMC simulations [^J. In the inset we report the correspond- 
ing data for £(T) derived from Eq. Q and from the spherical 
model as a function of the inverse temperature. According 
to the predictions of the QNLerM the QCR should occur for 
0.4 < T/J < 0.8 (i.e. 0.016<1/T<0.036), with a correspon- 
dent linear behavior of £(T) vs. 1/T, at variance with the 
experimental findings. 
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